Maximizing Spin Torque Diode Voltage by Optimizing Magnetization Alignment 
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The optimum condition of the magnetization alignment to maximize the spin torque diode voltage 
is derived by solving the Landau-Lifshitz-Gilbert equation. We show that the optimized diode 
voltage can be one order of magnitude larger than that of the conventional alignment where the 
easy axes of the free and the pinned layers are parallel. These analytical predictions are confirmed 
by numerical simulations. 



There has been great interest in spin-torque-induced 
magnetization dynamics [lHj] due to its potential appli- 
cation to spintronics devices such as magnetic random ac- 
cess memory (MRAM) and microwave oscillators. A spin 
torque diode [5l4l3|| is another important spintronics ap- 
plication, which enables us to rectify an alternating cur- 
rent in magnetic tunnel junction (MTJ) by synchroniz- 
ing the current with the resonant oscillation of the tunnel 
magnetoresistance (TMR). In 2010, a spin torque diode 
effect with relatively large sensitivity (~ 170 mV/mW) 
was observed experimentally ill]; however, the observed 
sensitivity was still lower than that of the Schottky diode. 

The spin torque diode effect arises from the combina- 
tion of the spin torque and the TMR effects. The spin 
torque originating from an alternating current induces 
a small oscillation of the magnetization of the free layer 
around its steady state, as a result of which the resistance 
of the MTJ oscillates through the TMR effect. The os- 
cillations of the current and the resistance create a direct 
voltage called a spin torque diode voltage. A large diode 
voltage, which determines the sensitivity of the diode, 
is obtained at the resonance frequency of the free layer. 
Hereafter, we refer to the MTJ in which the magnetiza- 



tion of the pinned layer is parallel to the easy axis of the 
free layer as the conventional alignment [5|, Iff . The diode 
voltage in the conventional alignment is on the order of 
10 - 100 fiV. A further increase of the diode voltage is 
desirable to excess the sensitivity of the Schottky diode. 

In this letter, we show that the spin torque diode volt- 
age can be significantly enhanced by choosing an appro- 
priate magnetization alignment of the free and pinned 
layers. We derive the optimum condition to maximize 
the diode voltage by solving the Landau-Lifshitz-Gilbert 
(LLG) equation. The optimum condition is determined 
by the competition between the contributions from the 
amplitude of the TMR oscillation and the linewidth of 
the power spectrum of the magnetization oscillation. We 
show that the optimum alignment shifts from the orthog- 
onal alignment, and that the diode voltage with the op- 
timized condition can be one order of magnitude larger 
than that of the conventional alignment. These analyti- 
cal predictions are confirmed by numerically solving the 
LLG equation. 

The system we consider is schematically shown in Fig. 
[TJ The MTJ consists of the free and the pinned layers 
separated by a nonmagnetic barrier. The x-axis is par- 



V 



(positive current) 



FIG. 1: A schematic view of the spin torque diode system. 
The positive current is defined as the electron flow from the 
free to the pinned layer. The easy axis of the free layer is 
parallel to the a;-axis. The unit vectors pointing in the direc- 
tion of the magnetizations of the free and the pinned layers 
are denoted as m and p, respectively. The angles tp, <pn, and 
ip-p are the steady state of the magnetization of the free layer 
and the directions of the applied field and the magnetization 
of the pinned layer, respectively. 



allel to the easy axis of the free layer while the z-axis is 
normal to the film plane. The unit vectors pointing in the 
direction of the magnetizations of the free and the pinned 
layers are denoted as m and p = (cos</?p,sin<y9 p ,0), re- 
spectively, where the direction of the magnetization of 
the pinned layer is described by the angle (p p . In the 
conventional alignment, tp p — 0. 

The magnetization dynamics of the free layer is de- 
scribed by using the LLG equation with the spin torque 
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Throughout this letter, 
magnetization dynamics 
the macrospin model. 
H = (ffappi cos ip H + Hunix 



we assume that 
is well-described 
The magnetic 
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i7 appl sin </?#,- 47tMto z ) is 
defined by the derivative the magnetic energy density, 
E = - M H app i(cos ip H m x + sin ip H m y ) - (MH K /2)ml + 
2irM 2 m 2 , , with respect to m, and consists of the applied 
field, -ffappi, the uniaxial anisotropy field along the 
easy axis, H^, and the demagnetization field along the 
hard axis, AttM. The gyromagnetic ratio and Gilbert 
damping constant are denoted as 7 and a, respectively. 
The spin torque p, |6| consists of the Slonczewski torque, 



a j, and the field like torque, bj, denned as 



aj 



2eMV ' 
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and bj = /3a j, respectively, where /3 is the ratio between 
the Slonczewski torque and the field like torque. Here, M 
and V are the saturation magnetization and the volume 
of the free layer, respectively. The spin polarization of 
the current is denoted as r\. The current / consists of 
the direct and the alternating currents, due to which aj 
and bj are decomposed into dc and ac parts as aj = 



flj(dc) + aj( ac ) 



and b 



J 



V(dc) 



'j(ac) 



respectively The 



positive current is defined as the electron flow from the 
free to the pinned layer. 

The spin torque diode effect arises from the small 
amplitude oscillation of the magnetization of the free 
layer around the steady state. Let us introduce 
two angles, (9,<p), characterizing the direction of 
the magnetization at the steady state as m^ ' = 
(sin 9 cos ip, sin 9 sin ip, cos 9). From eq. ([T|), the steady 
state should satisfy the following two conditions [17]: 

-Happi cos 9 cos((p — <ph) + (Hk cos 2 <p + AttM) sin 9 cos 9 
- aj(dc) sin(v3 - p p ) + 6 J(dc) cos9cos{ip - <p p ) = 0, 

(3) 



H, 



r a PP i sin(</3 -tp H ) + H K sin 9 sin tp cos ip 
- aj(dc) cos 9 cos((p - ip p ) - b J(dc) sm(<p - ip p ) = 0. 

In the absence of a direct current, the steady state is 
equal to the equilibrium state, i.e., the minimum state of 
the magnetic energy density E located in the film plane 
(9 = ir/2). Below, we set 9 = n/2 by assuming that the 
magnitudes of ajr dc \ and 6j(d c ) are small compared with 
|H| [lg|. Let us introduce the coordinates (X, Y,Z) in 
which the Z-axis is parallel to the steady state m^. The 
LLG equation can be linearized by applying the approx- 
imations mz — 1 and |mx|, \ttiy\ *C 1. The linearized 
LLG equation is given by 
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where we use the approximation 1 + a 2 ~ 1 because the 
Gilbert damping constant a is on the order of 10~ 2 [19j . 
The components of the coefficient matrix M are given by 

Mn = -a.i(dc)Pz + ol (H x + bj( dc )p z ) , (6) 



M12 = {H Y + b J{dc) pz) + aaj(dc)p z , (7) 



I21 = - (H x + b J(dc) p z ) - aa J{dc) p z , (8) 



I22 = -a,j(dc)Pz + a (H Y + bj {dc )p z ) , (9) 



where Hx and Hy are given by 

H x = -ffappi cos(y -p H )+H K cos 2 p + 4ttM, (10) 



Hy = iJ a ppi COs(i^ - <p H ) + H K COS 2lf. 



(11) 



The components of p in the XFZ-coordinate are 
given by {p x ,PY,Pz) = (0,-sin(<p - p p ),cos(<p> - ip p )). 
We assume that the alternating current is given by 
/ ac sin(27r/i) . Then, the solutions of mx and my can 
be obtained by solving eq. ([5]) . The explicit forms of mx 
and my , respectively, are given by 
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where 7 = 7/(271"), and aj( ac ) and 6j( ac ) are defined as 

aj(ac) = aj(ac) sin(27r/i) and 6j( ac) = 6,/( ac) sin(27r/t), re- 
spectively The resonant frequency / res and the linewidth 
A/ [2(| are defined as 

/res = 2~y (-Hx + &J(dc)Pz) (Hy + 6, /(dc) p Z ) + (a J(dc) p Z ) 2 , 

(14) 

A/ = ^ [a (JT X + flV) - 2a J{dc)P z] . (15) 

In the absence of a direct current, / ros is identical to the 
ferromagnetic resonant (FMR) frequency, /fmr- 
The spin torque diode voltage is defined as 
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I ac sin(27r/i), 
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where T = 1// and i? is the TMR given by R = 
i? P + (Aii/2)(1 - m • p), where Ai? = i? A p - i?p is the 
difference between the resistances at the parallel (Rp) 
and antiparallel (-Rap) alignments of the magnetizations. 
It should be noted that only the oscillation part of the 
TMR, — (Ai?/2)m ■ p, where mx and my oscillate with 
the frequency /, contributes to the diode voltage. By 
substituting eqs. (fT2)) and (fl3|) into eq. (fT6|) . the explicit 
form of Vdc is given by 
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FIG. 2: Dependence of the magnitude of the diode voltage, 
eq. (|20p . on the direction of the magnetization of the pinned 
layer. The values of the direct current are (a) 7dc = 0.5 mA 
and (b) I^ c — —0.5 mA. The conventional alignment corre- 
sponds to ipp = 0. The voltage is maximized at the optimized 
angle fp pt , and is zero at <p p = (p. 



Here, the Lorentzian and anti-Lorentzian parts, 
and srf(f), are given by 



#(J) = 



/ 2 A/7aj(ac)(l-p|) 
(/ 2 - fL) 2 + (/A/) 2 
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where Jf a = -aj( dc) pzPy and 3% = (H x +bj( dc) p z )p Y - 
The diode voltage, eq. ([17]). shows a peak near the 
resonant frequency / res . At / = / res , V dc is given by [1] 



Vdc(/res) 



AflJa 



aj(ac)Sin 2 (^-(^ p ) 



4 a(iTx + #r) - 2a J(dc ) cos(^ - f P ) 

where we use pz = cos(</? — ¥>p)- The term sin ((/? — <^ p ) 
in the numerator of eq. (|20[) arises from the oscilla- 
tion part of the TMR, —(AR/2)m ■ p, and is maxi- 
mized in the orthogonal alignment of the magnetizations, 
<p p — <p = 90° . The maximum diode voltage has been es- 
timated in this orthogonal alignment @. However, the 
diode voltage depends on not only the TMR but also 
the linewidth of the power spectrum of mx and my- 
The term 2aj/ dc \ cos(ip — (p p ) in the denominator of eq. 
(|2"U)) represents the enhancement or the reduction of the 
linewidth due to the spin torque acting as a damping 
or anti-damping factor, depending on the direction of 
the current. The optimum condition is determined by 
the competition between the contributions from the am- 
plitude of the TMR oscillation and the linewidth of the 
power spectrum of the magnetization oscillation. We find 
that the diode voltage, Vd c (/res), can be maximized when 
the magnetization of the pinned layer points to the di- 
rection 



<^° pt = cos" 1 




</>, 



(21) 



where the double sign "=p" means the upper (— ) for 
Idc/I c > and the lower (+) for 7<i c /J c < 0. The quan- 
tity I c is the absolute value of the critical current of the 
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FIG. 3: Spin torque diode voltages obtained by numerically 
solving the LLG equation, where the direction of the magneti- 
zation of the pinned layers and the value of the direct current 
are taken to be (a) (<p p ,/ d c(mA)) = (0,0.5), (b) (0,-0.5), (c) 
(<£° pt ,0.5), and (d) (v?° pt , -0.5), respectively. 

spin-torque-induced magnetization dynamics around the 
steady state defined as 



2aeMV 
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-ffappl COs((y9 - <fn) + H K 



cos 2 ip + cos 2ip 
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The optimum alignment, (p° pt 



( 22 ) 
.... ip, shifts from the or- 

thogonal alignment as long as the direct current is finite, 
while ip° pt —<p= 90° for I dc — because the spin torque 
does not affect the linewidth in this case. The maximized 
diode voltage, V^ c pt (/ res ), is given by 
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where the meaning of the double sign " =p " is the same 
as in eq. (121 [) . Equations ([21]) and (|23|) are the main 
results of this study. These results indicate that the spin 
torque diode voltage can be significantly enhanced by 
choosing an appropriate alignment of the magnetizations. 
It should be noted that since <p° pt is a real number, |ldc| 
should be less than |/ c |. In the opposite case, I dc /I c > 
1, the above formula is not applicable because the spin 
torque equals or overcomes the damping, due to which 
the steady state becomes unstable, and thus, the LLG 
equation cannot be linearized. In the limit of idc ~> 0, 
eq. d23) becomes ARI^J(8I C ). 

Let us quantitatively estimate how much the spin 
torque diode voltage can be enhanced by the optimiza- 
tion of the magnetization alignment. Figures EJa) and 
Hfb) show the dependence of the magnitudes of Vdc(/res) 
on the direction of the magnetization of the pinned layer, 
ip p , where Idc = 0.5 mA in (a) and —0.5 mA in (b). The 
values of the other parameters are taken to be M = 1000 
emu/c.c, H K = 200 Oe, ff appl = 100 Oe, tp H = 30°, 



y = 7rx80x35x2 nm 3 , 7 = 17.64 MHz/Oe, a = 0.01, 
7] = 0.5, p = 0.1, 7 ac = 0.2 mA, and AR = 100 O. The 
equilibrium direction of the free layer is estimated to be 
ip ~ 10° while the optimized directions of the pinned 
layer are <p° pt ~ 75° for I dc = 0.5 mA and 125° for 
7dc = —0.5 mA. It should be noted that the optimized 
alignment, (p° pt — ip, shifts from the orthogonal align- 
ment. In the conventional alignment (<p p = 0), the diode 
voltages are 72 //V for I dc = 0.5 mA and 13 /iV for 
Idc = —0.5 mA, where the diode voltage for I dc > 
is larger than that for I dc < because the positive cur- 
rent has an anti-damping effect in this case, and thus, 
reduces the linewidth. On the other hand, the magni- 
tude of the maximum diode voltage V^ is estimated to 
be 837 fiV for I dc = ±0.5 mA. Thus, the diode voltage 
satisfying the optimum condition is expected to be one 
order of magnitude larger than that of the conventional 
alignment. 

We confirmed the above analytical predictions by nu- 
merically solving the LLG equation [21j. Figures [3J a) 
and EJb) show the dependences of the diode voltages on 
the frequency of the alternating current at (p p = with 
(a) I dc = 0.5 and (b) —0.5 mA, respectively, while Figs. 
[1c) and^d) show V dc at <p p = <p° pt with (c) I dc = 0.5 
and (d) —0.5 mA. Although the direct current affects 
the resonant frequency, as shown in eq. (|14[) . the peak 
of the diode voltage appears approximately at the FMR 



frequency, /fmr = 5.3 GHz, because the magnitude of 
the direct current is relatively small; the peak frequen- 
cies of Fig. [3](a) and [3fb) are 5.4 GHz while those of 
Fig. [21[c) and[3][d) are 5.3 GHz. The magnitudes of the 
maximum voltage in Figs. [3ta) and [3jb) are 77 and 14 
/iV, while those of Figs. EJc) and E^d) are 900 and 897 
/iV, respectively. These results have a good agreement 
with Fig. [5J showing the validity of eqs. (|2"Tj) and (|23[) . 

In conclusion, we derived the optimum condition of the 
magnetization alignment of the free and the pinned layers 
to maximize the spin torque diode voltage by analyzing 
the competition between the oscillation of the tunnel- 
ing magnetoresistance and the reduction of the linewidth 
due to the spin torque. We showed that the optimum 
alignment shifts from the orthogonal alignment. We also 
showed that, under the optimized condition, the diode 
voltage can be one order of magnitude larger that that in 
the conventional alignment. These analytical predictions 
were confirmed by numerical simulations. The results 
indicate that the diode voltage can be significantly en- 
hanced by choosing an appropriate magnetization align- 
ment. Experimentally, the direction of the magnetization 
of the pinned layer maybe controlled during the anneal- 
ing process, as done in a TMR head [22J. 

The authors would like to acknowledge H. Kubota, H. 
Maehara, and S. Miwa for the valuable discussions they 
had with us. 



[1] J. C. Slonczewski, Phys. Rev. B 39, 6995 (1989). [13] 

[2] J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI 

(1996). 
[3] J. C. Slonczewski, J. Magn. Magn. Mater. 247, 324 [14] 

(2002). 
[4] L. Berger, Phys. Rev. B 54, 9353 (1996). [15] 

[5] A. A. Tulapurkar, Y. Suzuki, A. Fukushima, H. Kub- 
ota, H. Maehara, K. Tsunekawa, D. D. Djayaprawira, 

N. Watanabe, and S. Yuasa, Nature 438, 339 (2005). [16] 

[6] H. Kubota, A. Fukushima, K. Yakushiji, T. Naga- [17] 

hama, S. Yuasa, K. Ando, H. Maehara, Y. Nagamine, 

K. Tsunekawa, D. D. Djayaprawira, et al., Nature [18] 

Physics 4, 37 (2008). 
[7] J. C. Sankey, Y.-T. Cui, J. Z. Sun, J. C. Slonczewski, 

R. A. Buhrman, and D. C. Ralph, Nature Physics 4, 67 

(2008). 
[8] Y. Suzuki and H. Kubota, J. Phys. Soc. Jpn. 77, 031002 

(2008). 
[9] S. Yakata, H. Kubota, Y. Suzuki, K. Yakushiji, [19] 

A. Fukushima, S. Yuasa, and K. Ando, J. Appl. Phys. 

105, 07D131 (2009). 
[10] C. Wang, Y.-T. Cui, J. Z. Sun, J. A. Katine, R. A. [20] 

Buhrman, and D. C. Ralph, J. Appl. Phys. 106, 053905 

(2009). 
[11] S. Ishibashi, T. Seki, T. Nozaki, H. Kubota, S. Yakata, 

A. Fukushima, S. Yuasa, H. Maehara, K. Tsunekawa, [21] 

D. D. Djayaprawira, et al., Appl. Phys. Express 3, 073001 

(2010). [22] 

[12] S. Miwa, S.-Y. Park, S.-I. Kim, Y. Jo, N. Mizuochi, 

T. Shinjo, and Y. Suzuki, Appl. Phys. Express 5, 123001 

(2012). 



D. Bang, T. Taniguchi, H. Kubota, T. Yorozu, H. Ima- 
mura, K. Yakushiji, A. Fukushima, S. Yuasa, and 
K. Ando, J. Appl. Phys. Ill, 07C917 (2012). 

L. Landau and E. Lifshits, Phys. Z. Sowjetunion 8, 153 
(1935). 

E. M. Lifshitz and L. P. Pitaevskii, eds., Statistical 
Physics (Part 2), Course of Theoretical Physics, Vol. 9 
(Butterworth-Heinemann, Oxford, 1980), Chapter 7. 

T. L. Gilbert, IEEE Trans. Magn. 40, 3443 (2004). 
S. V. Vonsovskii, ed., Ferromagnetic Resonance (Perga- 
mon Press, Oxford, 1966), Chapter 2. 
In the experiments [fg|, the magnitude of the direct cur- 
rent is comparable to or smaller than that for the spin 
torque switching; see also eq. ((22}. It is well known that 
the critical current of the spin torque switching is on the 
order of a|H|. Since the Gilbert damping constant is very 
small (a ~ 10 -2 ) [l9|. we assume that the effect of the 
direct current on the determination of (0, ip) is negligible. 
M. Oogane, T. Wakitani, S. Yakata, R. Yilgin, Y. Ando, 
A. Sakuma, and T. Miyazaki, Jpn. J. Appl. Phys. 45, 
3889 (2006). 

We neglect the term fej(dc) in eq. (|15[1 because the field 
like torque is smaller than the Slonczewki torque, and 
the magnitude of the Slonczewki torque is on the order 
of a|H|, as mentioned in ref. [la ]. 

T. Taniguchi and H. Imamura, Appl. Phys. Express 4, 
103001 (2011). 

The magnetization alignment of a TMR head of a hard 
disk drive is set to be orthogonal by annealing process. 



